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Abstract. The mechanisms of information transmission are investigated in a 

lattice of coupled continuous maps, by analyzing the propagation of both finite 

£>v and infinitesimal disturbances. Two distinct regimes are detected: in the former 

' . case, both classes of perturbations spread with the same velocity; in the latter 

O 



case, finite perturbations propagate faster than infinitesimal ones. The transition 
between the two phases is also investigated by determining the scaling behaviour 
of the order parameter. 
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Recently, it has been discovered that coupled-map lattices (CML) can exhibit 
irregular behaviour even if the maximum Lyapunov exponent is negative [1, 
2]. This phenomenon has been called "stable chaos" [2] in order to stress the 
apparent incongruity. As a matter of fact, a periodic configuration is always 
approached in a finite time, on a lattice of length L. However, the average 
transient time T(L) diverges exponentially with L. In the thermodynamic limit 
L — ► oo, the statistical behaviour exhibited for times shorter than T(L) appears 
to be stationary. Therefore, the usual chaotic indicators (as, e.g., Lyapunov 
exponents) can be meaningfully applied to the characterization of the transient 
evolution as well. In analogy with standard transient chaos, the value of the 
Lyapunov exponents depends on the order the limits L — ► oo and t — ► oo are 
taken. 

The reference model considered in this Letter is a coupled-map lattice 
based on the usual diffusive-coupling scheme [3], 

xj+i = f{x\) (la) 

x\= £ -x\- 1 + {l-e)x\+ £ -x\ +1 (16) 

where t is the discrete time, i labels the lattice sites, e measures the coupling 
strength and periodic boundary conditions are assumed (x® = x\ , x t +1 = x\). 

In both [1] and [2] , the map / responsible for stable chaos is of the type 

{PiX + qi < x < x c 

l-(l-q 2 )(x-x c )/rj x c <x <x c + i] , (2) 

02 + Vl{x - X c - 7]) X c + 1] < X < 1 

(x c = (1 — qi)/pi) for rj — 0, i.e. it is discontinuous and piecewise linear. The 
single map considered in Ref. [1] (/i), corresponding to the parameter values 
Pi = P2 = 0.91, oi = .1, and 02 = 0, is characterized by a globally stable period- 
25 orbit. The parameters of the map in Ref. [2] (^2) are pi = 2.7, p 2 = 0.1, 
gi = and 02 = 0.07. In such a case, a globally stable period-3 orbit exists. 

One can easily conjecture that the sudden amplification of the distance 
between points lying across the discontinuity (a phenomenon which escapes the 
linear stability analysis) is at the origin of the random evolution. In fact, Ershov 



and Potapov [4] have shown that nonchaotic turbulence persists in the CML if 
the map f\ is made continuous by introducing a finite r\ < 4 • 10 -4 . However, 
they have also shown that the continuization of the map destabilizes the period- 
25 orbit, leading to a standard chaotic attractor in the single map for rj > 
5 ■ 10 -3 Therefore, they have concluded that the observed irregular behaviour is 
a reminiscence of the usual transient chaos, characterized by a positive Lyapunov 
exponent. 

We have found that, in analogy to what happens for the model of Ref. [1], 
the dynamics of the CML with a continuous ji exhibits a negative maximum 
Lyapunov exponent if rj < r\\ — 4 • 10 -4 . However, in our case, the critical 
value of rj corresponding to the destabilization of the period-3 orbit is much 
higher (772 = P1Q2 — 1/pi — 0.14). Therefore, the phenomenon of stable chaos 
has a less obvious origin than conjectured in Ref. [4]. Nevertheless, one must 
conclude that in both models, the linear stability analysis is not able to capture 
the essence of the observed irregular behaviour. 

We conjecture that two different mechanisms of information flow on the 
lattice must be invoked: one based on the usual (local) linear instability, the 
other based on "nonlinear" propagation phenomena, analogous to those present 
in cellular automata [5] . The aim of the present Letter is to perform a compar- 
ative study of the two mechanisms. 

The CML has been analyzed for e = 2/3. We investigate the evolution of 
an initially localized disturbance u^ =0 = 5^o in the linear approximation. Its 
time evolution is of the type 
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3 A(i/t)t ^ (3) 



where the comoving Lyapunov exponent A(v) [6] is the growth rate of the distur- 
bance in a reference frame moving with velocity v = i/t. Within the light-cone 
defined by 

A(u £ ) = , 

the perturbation is exponentially amplified. Accordingly, vg represents the limit 
velocity of propagation of disturbances in the linear approximation. The direct 



estimate of A(v) may be affected by two technical problems: (i) a site-dependent 
normalization of the amplitude \w\\ to cope with the different growth rates along 
the different lines i = vt; (ii) the increasing number of sites that must be up- 
dated. These problems can lead to serious difficulties, whenever the convergence 
towards the asymptotic shape is slow. Deissler and Kaneko [6] applied the usual 
algorithm for Lyapunov exponents in a moving reference frame. This approach 
is not affected by the above problems and it works perfectly in open-flow sys- 
tems. However, in closed systems (like those ones we are currently investigating) 
boundary conditions strongly affect the evolution of Lyapunov vectors. 

This further problem can be overcome by following the approach recently 
developed in [7]. The authors have introduced the specific Lyapunov exponent 
A(/x), defined as the maximum growth rate of a perturbation w\ characterized 
by an exponential profile, w\ = exp (— fxi) ■ u\. The new variable u\, assumed to 
satisfy periodic boundary conditions, obeys the evolution equation 

< +1 = f'(yl)[^ur 1 + (l-e)ul + £ -e^ui +1 } • (4) 

The specific exponent A(/i) can be estimated by applying the standard algorithm 
developed for the computation of Lyapunov spectra. The comoving exponent is 
related to A(//) through the Legendre trasform [7] 

A(v) = A(» - /j\'(/j) ; v = -A'(» (5) 

From the above equation, 

v e = -A' (71) , 

where fl is the /(/-value such that the first of Eq. (5) vanishes. The derivative 
A'(/i) can be estimated with a good accuracy by iterating a suitable recursive 
equation. After deriving Eq. (4) with respect to the parameter \i and introducing 

z\ = (u l t )', we obtain 

4+i = / / (yD[-|e-^r 1 +|e^^ 1 ]+/ / (yi)[|e-^r 1 +(l- £ )^+|e^|+ 1 ] (6) 
The derivative of the specific exponent is then given by 

W (* + i)(£ fe ^ +1 -^ + i) 1/2 {) 
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The results of simulations performed for several rj values are reported in Fig. 1 
(dashed line). The linear propagation develops for rj > rj\, when the maximum 
Lyapunov exponent becomes larger than zero. 

The evolution of finite disturbances and the related nonlinear mechanisms 
can be investigated only in terms of a direct approach. Let us consider two 
initial configurations {x l t }, {y l t } differing by order 0(1) in a finite interval and 
coinciding outside. The corresponding configurations, generated after T time 
steps, are then compared to determine the left and right borders ii, resp. i r of the 
region where they differ more than a preassigned threshold 9. The propagation 
velocity v n g of finite disturbances is then defined as 

,. (i r -ii) 
v ne = hm 



Numerical simulations reveal a sizeable dependence of the actual value of v n £ on 
the computational accuracy. For instance, upon increasing the number of signif- 
icant digits from 8 to 64, the value of v n £ increases by about 2 % for rj = 0.01. 
This inaccuracy prevents a quantitative study of the scaling behaviour of the 
difference Av — V£ — v n £ for V£ — ► v n £. Such a difficulty has been circumvented 
by extrapolating the numerically exact value of v n £ from the results of increas- 
ingly accurate simulations [8]. As a matter of fact, the finite-precision estimates 
of the velocity converge as 1/N%, where N\, is the number of significant digits. 
The asymptotic values of v n £ are reported in Fig. 1 (solid line). Moreover, upon 
varying the value of the threshold 9 from 1.10 -4 to 0.1, the resulting variation 
of v n £ is practically negligible. 

Any infinitesimal perturbation w\ is exponentially amplified within the 
light-cone defined by V£. Therefore, w\ eventually becomes larger than any 
finite threshold 9, so that the nonlinear velocity v n g cannot be smaller than 
V£. This inequality is indeed satisfied in the simulations reported in Fig. 1, 
where two different "phases" are identified: for rj > rj c = 0.013, V£ = v n £ (I); 
for rj < Tj c , v n £ > V£ (II). The former phase corresponds to the usual regime 
found in lattices of logistic and cubic maps [9] . The latter phase corresponds to 
a new regime, where the instability due to a positive Lyapunov exponent is not 
sufficient to account for the propagation of information along the chain. The 



examples analysed in Refs. [1,2] represent a limit case of this scenario, where 
there is no linear instability at all. 

The transition between the two regimes can be investigated from the de- 
pendence of the order parameter Av = vg — v n g versus 1] C — rj. The numerical 
results reported in Fig. 2, reveal a scaling behaviour of the type 

Av = (r] c - r,y . (8) 

The optimal scaling behaviour is obtained for the above reported value of r\ c = 
0.013 yielding 7 ~ 2. Having we failed to derive even a heuristic explanation 
of this exponent, we have tried to reach a more complete understanding of this 
phemomenon by studying the fluctuations of the effective maximal comoving 
Lyapunov exponent A. The effective exponent A corresponding to a velocity v 
over a time t is defined as 



A = - lot 
t 



Vt ~ x t 



(9) 



where i = vt and {x\}, {y l t } represent two configurations which, at time t = 0, 
differ by Sq = 5 in the site i — 0, while coincide elsewhere. Moreover, let 
P(A,v,t) denote the probability density to find A in the interval (A, A + dA). 
In the limit 5^0 and t — ► 00, the whole information about the fluctuations of 
the Lyapunov exponent is conveyed by the multifractal spectrum 

s(kv) = ^(nkvA) (10) 

I.' 

The Lyapunov analysis has been performed for 77 = 2 • 10 -3 . For this parameter 
value, linear and nonlinear velocities are vg = 0.4184, v n g = 0.5805, respectively. 
We have computed the spectrum S for different velocities and time lags. The 
curves reported in Fig. 3 have been all obtained for t = 40; labels 1, 2 and 3 
refer to 5 = 10 -7 , 10 -4 and 10 -2 , respectively. Finally, the spectra reported in 
Fig. 3a correspond to v — 0.7, while those ones in Fig. 3b to v = 0.5. 

If the amplitude of the perturbation is initially smaller than S c = 
exp[— A max (v)t] (where A max (v) is the maximum value assumed by the effective 
Lyapunov exponent A at velocity v), then it will always obey a linear equation 
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during its evolution over the first (arbitrary) t time steps. In this case, one is 
simply estimating the multifractal spectrum of the comoving Lyapunov exponent 
A. On very general grounds, a bell shaped spectrum is expected. This is in fact 
what observed for the curves labelled with "1" in Fig. 3. If 5 is increased above 
S c , some perturbations become of order 0(1) and are then controlled by the full 
nonlinear equations. One could simply conjecture that all A values greater than 
A u = — log S/t are no longer detectable, while the rest of the spectrum remains 
unchanged. This is precisely what happens when the velocity v is larger than 
v n £ as, for instance, in Fig. 3a. If instead, v < v n £, a second peak appears at high 
A- values, when S is increased (see Fig. 3b). For 6 sufficiently "large", the new 
peak becomes higher than the old one, meaning that the nonlinear mechanism 
prevails onto the linear one. 

The order of magnitude of the statistical error is much less than the ampli- 
tude of the oscillations observed at small As: they are intrinsic fluctuations which 
are seen to decay slowly for increasing t. However, this finite-size phenomenon 
does not prevent a clear-cut distinction between the qualitative behaviour of 
the spectra observed in Fig. 3a and 3b. Therefore, we can conclude that the 
exchange of the two limits 5^0 and t — ► oo plays a crucial role in the phe- 
nomenon studied in this Letter. The usual definition of Lyapunov exponent is 
recovered if the former limit is taken first, while the nonlinear propagation of 
information is revealed only by taking first the latter limit. 

In this Letter we have discussed two different mechanisms for information 
propagation. We have identified a regime where the prevailing mechanisms de- 
pends on the amplitude of the propagating perturbation. Whenever nonlinear 
processes eventually overtake the linear amplification, we expect different statis- 
tical properties of the associated spatio-temporal pattern. Finally, notice that 
the relevance of nonlinear mechanisms has been pointed out also in a standard 
CML of logistic maps, with reference to the predictability problem [10]. 
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Figure Captions 

Fig. 1 Linear (dashed) and nonlinear (solid) velocities versus rj. For the CML 
model (2) with p\ = 2.7,£>2 = 0.1, q± = , qi = 0.07 and diffusive coupling 
e = 2/3. The two arrows indicate the position of the parameters r/i = 
4 • 10" 4 and rj c = 0.013. 

Fig. 2 Av = V£ — v n £ versus r\ c — rj (rj c = 0.013) reported in a log - log scale. 
The circles represent the actual point where the difference Av have been 
calculated. 

Fig. 3 Multifractal spectra S(A,v) for rj = 2 • 10" 3 and t = 40: (a) v = 0.7; (b) 
v = 0.5. The labels 1, 2 and 3 correspond to 5 = 10~ 7 , 10~ 4 and 10~ 2 , 
respectively. The spectra have been obtained as an average over 200,000 
different realizations. 



